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In this final lecture, we will study the behaviour of the Bloch equations in different regimes of 
resonance and relaxation. The Bloch equations are formulated as a vector model, and numerical 
solutions to the equations are discussed. 

For steady-state interaction, the polarization density of the medium, as obtained from the 
Bloch equations, is expressed in a closed form. The closed solution is then expanded in a power 
series, which when compared with the series obtained from the susceptibility formalism finally tie 
together the Bloch theory with the susceptibilities. 

The outline for this lecture is: 

• Recapitulation of the Bloch equations 

• The vector model of the Bloch equations 

• Special cases and examples 

• Steady-state regime 

• The intensity dependent refractive index at steady-state 

• Comparison with the susceptibility model 

Recapitulation of the Bloch equations for two-level systems 

Assuming two states \a) and \b) to be sufficiently similar in order for their respective lifetimes T a ss 
T b « T\ to hold, where T\ is the longitudinal relaxation time, the Bloch equations for the two-level 
are given as 

— =-Av-u/T 2 , (la) 
dv 

— =Au + f3(t)w~v/T 2 , (lb) 
^ = -0(t)v-(w-w o )/T 1 , (lc) 

where (5 = er^ b E"(t)/h is the Rabi frequency, being a quantity linear in the applied electric field 
of the light, A = Q ba — w is the detuning of the angular frequency of the light from the transition 
frequency Q ba = (E b — K a )/h, and where the variables (u 7 v 7 w) are related to the matrix elements 
p mn of the density operator as 

u = Pba + Pabi 

v = l {Pba- Pab), 
W — Pbb - Paa- 

In these equations, p^ b is the temporal envelope of the off-diagonal elements, given by 

Pab = Pa b exp[i(fl ba - A)t]. 

In the Bloch equations (1), the variable w describes the population inversion of the two-level 
system, while u and v are related to the dispersive and absorptive components of the polarization 
density of the medium. In the Bloch equations, wo = po(b) — Po(a) is the thermal equilibrium 
inversion of the system with no optical field applied. 
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The resulting electric polarization density of the medium 

The so far developed theory of the density matrix under resonant interaction can now be applied 
to the calculation of the electric polarization density of the medium, consisting of N identical 
molecules per unit volume, as 

P^v,t) = N{ef^) 

= NTr[per^] 

k—a.b 

= N H E ( k \P\j)(j\ef,\k) 

k—a.b j—a,b 

= N J2 {(k\p\a){a\ef„\k) + {k\p\b)(b\er^\k)} 

k—a.b 

= N{(a\p\a) <a|ef» + (b\p\a)(a\ef^\b) + (a\p\b)(b\ef^\a) + (b\p\b)(b\ef^\b)} 
= N(p ba er% b + Paber£ a ) 

— {Make use of p ab = (u + iv) exp(iuit) — p* ha ) 

= N[(u — iv) exp(—iLvt)er^ b + (u + iv) exp(iwi)er^ a ]. 

The temporal envelope P£ of the polarization density is throughout this course as well as in Butcher 
and Cotter's book taken as 

P"(r,t) = Re[P^exp(-twt)] ) 

and by identifying this expression with the right-hand side of the result above, we hence finally 
have obtained the polarization density in terms of the Bloch parameters (u,v,w) as 

P£{T,t) = Ner2 b (u-iv). (2) 

This expression for the temporal envelope of the polarization density is exactly in the same mode 
of description as the one as previously used in the susceptibility theory, as in the wave equations 
developed in lecture eight. The only difference is that now we instead consider the polarization 
density as given by a non-perturbative analysis. Taken together with the Maxwell's equations (or 
the propcrn wave equation for the envelopes of the fields), the Bloch equations are known as the 
Maxwell-Bloch equations. 

From Eq. (2), it should now be clear that the Bloch variable u essentially gives the in-phase 
part of the polarization density (at least in this case, where we may consider the transition dipolc 
moments to be real- valued) , corresponding to the dispersive components of the interaction between 
light and matter, while the Bloch variable v on the other hand gives terms which are shifted ninety 
degrees out of phase with the optical field, hence corresponding to absorptive terms. 
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The vector model of the Bloch equations 

In the form of Eqs. (1), the Bloch equations can be expressed in the form of an Euler equation as 

— =nxTL-(u/T 2 ,v/T 2 ,(w-w )/T 1 ), [B.&C. (6.54)] 

at v ^ / 

relaxation term 

where R = (u, v, w) is the so-called Bloch vector, that in the abstract (e„, e„, e w )-space describes 
the state of the medium, and 

n = (-/?(t),o,A) 

is the vector that gives the precession of the Bloch vector (see Fig. 1). 

This form, originally proposed in 1946 by Felix Bloch 1 for the motion of a nuclear spin in 
a magnetic field under influence of radio-frequency electromagnetic fields, and later on adopted 
by Feynman, Vernon, and Hellwarth 2 for solving problems in maser theory 3 , corresponds to the 
motion of a damped gyroscope in the presence of a gravitational field. In this analogy, the vector 
fi can be considered as the torque vector of the spinning top of the gyroscope. 

w 

A 



v 



Roc = (0,0, w ) 

Figure 1. Evolution of the Bloch vector R(t) = (u(t),v(t),w(t)) around the "torque vector" 
fi = (—/3(f), 0, A). In the absence of optical fields, the Bloch vector relax towards the thermal 
equilibrium state Roo = (0, 0,Wq), where wq — p(b) — p(a) is the molecular population inversion 
at thermal equilibrium. At moderate temperatures, the thermal equilibrium population inversion 
is very close to w = —1. 

From the vector form of the Bloch equations, it is found that the Bloch vector rotates around 
the torque vector 17 as the state of matter approaches steady state. For an adiabatically changing 
applied optical field (i. e. a slowly varying envelope of the field), this precession follows the torque 
vector. 

The relaxation term in the vector Bloch equations also tells us that the relaxation along the 
w-direction is given by the time constant T l7 while the relaxation in the (u, w)-plane instead is 
given by the time constant T 2 . By considering the w-axis as the "longitudinal" direction and the 
(u, w)-plane as the "transverse" plane, the terminology for T\ as being the "longitudinal relaxation 
time" and T 2 as being the "transverse relaxation time" should hence be clear. 

1 F. Bloch, Nuclear induction, Phys. Rev. 70, 460 (1946). Felix Bloch was in 1952 awarded the 
Nobel prize in physics, together with Edward Mills Purcell, "for their development of new methods 
for nuclear magnetic precision measurements and discoveries in connection therewith" . 

2 R. P. Feynman, F. L. Vernon, and R. W. Hellwarth, Geometrical representation of the Schrdd- 
inger equation for solving maser problems, J. Appl. Phys. 28, 49 (1957). 

3 Microwave Amplification by Stimulated Emission of Radiation, a device for amplification of 
microwaves, essentially working on the same principle as the laser. 
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Transient build-up at exact resonance as the optical field is switched on 

The case T\ » T 2 - Longitudinal relaxation slower than transverse relaxation 
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Figure 2a. Evolution of the Bloch vector (u(t), v(t), w(t)) as the optical field is switched on, 
for the exactly resonant case (S = 0), and with the longitudinal relaxation time being much 
greater than the transverse relaxation time (T\ 3> T2). The parameters used in the simulation are 
r) = Ti/T 2 = 100, S = AT 2 = 0, w = -1, and j(t) = /3(t)T 2 = 3, t > 0. The medium was initially 
at thermal equilibrium, (u(0),v(0),w(0)) = (0,0, wq) — —(0,0,1). 




x= tl T 2 

Figure 2b. Evolution of the magnitude of the polarization density |P w (t)| ~ \u(t) — iv(t)\ as the 
optical field is switched on, corresponding to the simulation shown in Fig. 2a. 
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The case T\ « T 2 - Longitudinal relaxation approximately equal to transverse relaxation 
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Figure 3a. Evolution of the Bloch vector (u(t),v(t),w(t)) as the optical field is switched on, for 
the exactly resonant case (5 = 0), and with the longitudinal relaxation time being approximately 
equal to the transverse relaxation time (T\ « T 2 ). The parameters used in the simulation are 
r) = Ti/T 2 = 2,6 = AT 2 =0,w = -1, and j(t) = /3(t)T 2 = 3, t > 0. The medium was initially 
at thermal equilibrium, (u(0),v(0),w(0)) = (0,0, wq) = —(0,0,1). 




x= tl T 2 

Figure 3b. Evolution of the magnitude of the polarization density |i"L(i)| ~ \u(t) — iv(t)\ as the 
optical field is switched on, corresponding to the simulation shown in Fig. 3a. 
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Transient build-up at off-resonance as the optical field is switched on 

The case T\ « T 2 - Longitudinal relaxation approximately equal to transverse relaxation 
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Figure 4a. Evolution of the Bloch vector (u(t), v(t), w(t)) as the optical field is switched on, 
for the off-resonant case (6^0), and with the longitudinal relaxation time being approximately 
equal to the transverse relaxation time (T\ « T 2 ). The parameters used in the simulation are 
n = Ti/T 2 = 2,6 = AT 2 = 4, w = -1, and j(t) = 0(t)T 2 = 3, t > 0. The medium was initially 
at thermal equilibrium, (u(0),v(0),w(0)) = (0,0, wq) = —(0,0,1). 




x= tl T 



Figure 4b. Evolution of the magnitude of the polarization density |P w (t)| ~ \u(t) — iv(t)\ as the 
optical field is switched on, corresponding to the simulation shown in Fig. 4a. 
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Transient decay for a process tuned to exact resonance 

The case Ti » T 2 - Longitudinal relaxation slower than transverse relaxation 
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Figure 5. Evolution of the Bloch vector (u(t),v(t),w(t)) after the optical field is switched off, for 
the case of tuning to exact resonance ((5 = 0), and with the longitudinal relaxation time being much 
greater than the transverse relaxation time (Ti 3> T 2 ). The parameters used in the simulation arc 
7/ = Ti/T 2 = 100, 5 = AT 2 = 0, w = -1, and 7 (i) = (3{t)T 2 = 0. 

The case T\ « T 2 Longitudinal relaxation approximately equal to transverse relaxation 



S o 



-0.5 - 




Figure 6. Evolution of the Bloch vector (u(t),v(t),w(t)) after the optical field is switched off, for 
the case of tuning to exact resonance (5 — 0), and with the longitudinal relaxation time being 
approximately equal to the transverse relaxation time (Ti ~ T 2 ). The parameters used in the 
simulation arc 77 = T 1 /T 2 = 2,6 = AT 2 = 0, w a = -1, and j(t) = (3(t)T 2 = 0. 
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Figure 7a. Evolution of the Bloch vector (u(t),v(t),w(t)) after the optical field is switched off, for 
the off- resonant case (6 ^ 0) , and with the longitudinal relaxation time being much greater than the 
transverse relaxation time (T\ 3> T2). The parameters used in the simulation are rj = T1/T2 = 100, 
S = AT 2 = 2, w = — 1, and j(t) = /3{t)T 2 — 0. (Compare with Fig. 5 for the exactly resonant 
case.) 



8 10 12 14 16 18 



8 10 12 14 16 18 

tlT 2 



Figure 7b. Same as Fig. 7a, but with 5 = — 2 as negative. 
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The case T\ « T 2 - Longitudinal relaxation approximately equal to transverse relaxation 
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Figure 8a. Evolution of the Bloch vector (u(t), v(t), w(t)) after the optical field is switched off, 
for the off-resonant case (5^0), and with the longitudinal relaxation time being approximately 
equal to the transverse relaxation time (T\ « T 2 ). The parameters used in the simulation are 
r) = Ti/T 2 = 2, 5 = AT 2 = 2, w Q = -1, and j(t) = /3(t)T 2 = 0. (Compare with Fig. 6 for the 
exactly resonant case.) 




x= tl T 2 

Figure 8b. Evolution of the magnitude of the polarization density ~ \u(t) — iv(t)\ as the 

optical field is switched on, corresponding to the simulation shown in Fig. 8a. 
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Figure 9a. Evolution of the Bloch vector (u(t),v(t),w(t)) after the optical field is switched off, for 
the far off- resonant case (5 ^ 0), and with the longitudinal relaxation time being much greater 
than the transverse relaxation time (Ti ;§> T 2 ). The parameters used in the simulation are 77 = 
T1/T2 = 100, 5 = AT 2 = 20, w Q = -1, and 7 (i) = f3(t)T 2 = 0. (Compare with Fig. 5 for the 
exactly resonant case, and with Fig. 7a for the slightly off-resonant case.) 




x= tl T 2 

Figure 9b. Evolution of the magnitude of the polarization density ~ \u(t) — iv(t)\ as the 

optical field is switched on, corresponding to the simulation shown in Fig. 9a. 
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The case T\ « T 2 - Longitudinal relaxation approximately equal to transverse relaxation 
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Figure 10a. Evolution of the Bloch vector (u(t), v(t), w(t)) after the optical field is switched off, for 
the far off-resonant case (6^0), and with the longitudinal relaxation time being approximately 
equal to the transverse relaxation time (Ti « T 2 ). The parameters used in the simulation are 
r) = Tx/T 2 =2,6 = AT 2 = 20, w = -1, and ~/(t) = /3(t)T 2 = 0. (Compare with Fig. 6 for the 
exactly resonant case, and with Fig. 8a for the slightly off-resonant case.) 




x= tl T 2 

Figure 10b. Evolution of the magnitude of the polarization density |-FL(£)| ~ — iv(t)\ as the 
optical field is switched on, corresponding to the simulation shown in Fig. 10a. 
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The case T\ <C T 2 - Longitudinal relaxation faster than transverse relaxation 
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Figure 11a. Same parameter values as in Fig. 6, but with the longitudinal relaxation time being 
much smaller than the transverse relaxation time (Ti <C T 2 ), n = Ti/T 2 = 0.1. (Compare with 
Figs. 9a and 10a for the cases T x T 2 and T\ ~T 2 , respectively.) 




x= tl T 
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Figure lib. Evolution of the magnitude of the polarization density |-P w (i)| ~ \u(t) — iv(t)\ as the 
optical field is switched on, corresponding to the simulation shown in Fig. 11a. 
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The connection between the Bloch equations and the susceptibility 

As an example of the connection between the polarization density obtained from the Bloch equa- 
tions and the one obtained from the susceptibility formalism, we will now - once again - consider 
the intensity-dependent refractive of the medium. 

The intensity- dependent refractive index in the susceptibility formalism 

Previously in this course, the intensity-dependent refractive index has been obtained from the 
optical Kerr-effect in isotropic media, in the form 

n = n + n 2 |E a) | 2 , 
where n = [1 + \xx (— uj; uj)] 1 / 2 is the linear refractive index, and 

is the parameter of the intensity dependent contribution. However, since we by now are fully aware 
that the polarization density in the description of the susceptibility formalism originally is given 
as an infinity series expansion, we may expect that the general form of the intensity dependent 
refractive index rather would be as a power series in the intensity, 

n = n a + n 2 |E w | 2 + n 4 |E w | 4 + n 6 |E w | 6 + . . . 

For linearly polarized light, say along the a;- axis of a Cartesian coordinate system, we know that 
such a series is readily possible to derive in terms of the susceptibility formalism, with the different 
order terms of the refractive index expansion given by the elements 

n 6 ~ Xxxxxxxxx(-U\ V, U, -UJ, UJ, -UJ, UJ, -Uj), 



Such an analysis would, however, be extremely cumbersome when it comes to the analysis of 
higher-order effects, and the obtained sum of various order terms would also be almost impossible 
to obtain a closed expression for. For future reference, to be used in the interpretation of the 
polarization density given by the Bloch equations, the intensity dependent polarization density is 
though shown in its explicit form below, including up to the seventh order interaction term in the 
Butcher and Cotter convention, 

K = eoX$ (-w; ")K (order n = 1) 

+ e (S/4) X i 3 l x (-uj;uj,uj,-uj)\EZ\ 2 EZ (order n = 3) 

+ e (5/8) X ( xL xxx (-uj; uj,uj, -uj,uj, -uj)\EZ\ 4 EZ (order n = 5) 

+ e {35/64:)x x % xxxxx (-uj;uj,uj, -u,u, -uj,uj, -u))\E^\ 6 E^ ( order n = 1 ) 
+ ... 

The other approach to calculation of the polarization density, as we next will outline, is to use the 
steady-state solutions to the Bloch equations. 
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The intensity-dependent refractive index in the Bloch-vector formalism 

For steady-state interaction between light and matter, the solutions to the Bloch equations yield 

u ~ iv = A-^/TV [B.&C. (6.53a)], 

Wo [l + (AT 2 ) 2 ] 

w=- ,. m „ 7^ ~ , B.&C. (6.53b) , 

where, as previously, (3 = er" h E"(t)/h is the Rabi frequency, though now considered to be a 
slowly varying (adiabatically following) quantity, due to the assumption of steady-state behaviour. 
From the steady-state solutions, the ^-component (/i = x, y, z) of the electric polarization density 
P(r, t) = Re[P w exp(— iuit)] of the medium hence is given as 

PS = Ner^(u-iv) 

= - Ner » ? w 

ab A-i/T 2 

= _ Ner » P w [l + {AT 2 f] (3) 

ab (A - i/T 2 ) [1 + (AT 2 )2 + /32 Tl T 2 ] 

r M ft 
= -New ^ 



(A - i/T 2 ) 



i i Zili A2 

1 + (1 + (AT 2 )2)P 



In this expression for the polarization density, it might at a first glance seem as it is negative for 
a positive Rabi frequency [3, henc giving a polarization density that is directed anti-parallel to 
the electric field. However, the quantity w = po(b) — po( a ), the population inversion at thermal 
equilibrium, is always negative (since we for sure do not have any population inversion at thermal 
equilibrium, for which we rather expect the molecules to occupy the lower state), hence ensuring 
that the off-resonant, real- valued polarization density always is directed along the direction of the 
electric field of the light. 

Next observation is that the polarization density no longer is expressed as a power series in 
terms of the electric field, but rather as a rational function, 

P^X/(1 + X*), (4) 

where 

X = VTiT 2 /(l + (AT 2 )2)/3 

= VTiT 2 /(l + (AT 2 )2)erM(t)/fi 

is a parameter linear in the electric field. The principal shape of the rational function in Eq. (4) is 
shown in Fig. 12. 

From Eq. (4), the polarization density is found to increase with increasing X up to X = 1, as 
we expect for an increasing power of an optical beam. However, for X > 1, we find the somewhat 
surprising fact that the polarization density instead decrease with an increasing intensity; this 
peculiar suggested behaviour should hence be explained before continuing. 

The first observation we may do is that the linear polarizability (i. e. what we usually associate 
with linear optics) follows the first order approximation p(X) = X. In the region where the 
peculiar decrease of the polarization density appear, the difference between the suggested nonlinear 
polarization density and the one given by the linear approximation is huge, and since we a priori 
expect nonlinear contributions to be small compared to the alsways present linear ones, this is 
already an indication of that we in all practical situations do not have to consider the descrease of 
polarization density as shown in Fig. 12. 

For optical fields of the strength that would give rise to nonlinearities exceeding the linear 
terms, the underlying physics will rather belong to the field of plasma and high-energy physics, 
rather than a bound-charge description of gases and solids. This implies that the validity of the 
models here applied (bound charges, Hamiltonians being linear in the optical field, etc.) are limited 
to a range well within X < 1. 



14 



Lecture notes 12 



Nonlinear Optics 5A5513 (2003) 



i 

i 

i 




X 

Figure 12. The principal shape of the electric polarization density of the medium, as function of 
the applied electric field of the light. In this figure, X = yjT x T 2 l{\ + (AT 2 ) 2 )0 is a normalized 
parameter describing the field strength of the electric field of the light. 



Another interesting point we may observe is a more mathematically related one. In Fig. 12, we 
see that even for very high order terms (such as the approximating power series of degree 31, as 
shown in the figure), all power series expansions fail before reaching X = 1. The reason for this is 
that the power series that approximate the rational function X/(l + X 2 ), 

X/(l + X 2 ) = X - X 3 + X b - X 7 + 

is convergent only for \X\ < 1; for all other values, the series is divergent. This means that no 
matter how many terms we include in the power series in X, it will nevertheless fail when it comes 
to the evaluation for \X\ > 1. Since this power series expansion is equivalent to the expansion 
of the nonlinear polarization density in terms of the electrical field of the light (keeping in mind 
that X here actually is linear in the electric field and hence strictly can be considered as the field 
variable) , this also is an indication that at this point the whole susceptibility formalism fail to give 
a proper description at this working point. 

This is an excellent illustration of the downturn of the susceptibility description of interaction 
between light and matter; no matter how many terms we may include in the power series of the 
electrcal field, it will at some point nevertheless fail to give the total picture of the interaction, and 
we must then instead seek other tools. 

Returning to the polarization density given by Eq. (3), we may now express this in an explicit 
form by inserting A = Q b a — v I0r the angular frequency detuning, the Rabi frequency = 
er° b E"(t)/h, and the thermal equilibrium inversion w = po(b) — p (a). This gives the polarization 
density of the medium as 



Np 2 r M r a 

PZ = So—rr(p (a)-p (b)Y 



e h (Q ba -u> - i/T 2 ) 



Till ( pr 1 E?T/fc\S 



x\±<L{~ f° r a two level medium 



nonlinear correction factor to Xmc*(~ uj '^ lj ) 



the field corrected susceptibility, xt^jE^) [sec Butcher and Cotter, section 6.3.1] 
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In this form, the polarization density is given as the product with a term which is identical to the 
linear susceptibility 4 (as obtained in the perturbation analysis in the frame of the susceptibility 
formalism), and a correction factor which is a nonlinear function of the electric field. 

The nonlinear correction factor, of the form 1/(1 + X 2 ), with X = \jT\T 2 j(}- + (AT 2 ) 2 )f3 as 
previously, can now be expanded in a power series around the small-signal limit X = 0, using 



1/(1 + X 1 ) = 1 - X z + X 4 - x b + . . . , 

from which we obtain the polarization density as a power series in the electric field (which for the 
sake of simplicitly now is taken as linearly polarized along the x-axis) as 

iVe 4 r x r x (r x ) 2 2 

soh^ (Q ba -w - 1/T2) [1/T.f + (Q ba - w) 2 J(T 2 /Ti) (5) 



iVe 6 r x r x (r x ) 4 4 

e h b (iiba -u- i/T 2 ) [i/Tj + (ii ba - uy\ z {T 2 /Tiy 

+ ... 

This form is identical to one as obtained in the susceptibility formalism; however, the steps that 
led us to this expression for the polarization density do not rely on the perturbation theory of 
the density operator, but rather on the explicit form of the steady-state solutions to the Bloch 
equations. 

Summary of the Bloch and susceptibility polarization densities 

To summarize this last lecture on the Bloch equations, expressing the involved parameters in the 
same style as previously used in the description of the susceptibility formalism, the polarization 
density obtained from the steady-state solutions to the Bloch equations is 

ive ' ' ^ 1 ab' ab ' 



P% = e — { Po {a) - p (b)) 



e h (Q ba - u - i/T 2 ) 



1 + MCT?)( er ») 2 



E 



By expanding this in a power series in the electrical field, one obtains the form (5), in which we 
from the same description of the polarization density in the susceptibility formalism can identify 

x W(- u . u ) = —(p (a)-p (b))-. 



e h (Qba - w - i/T 2 y 

AMpi r x r x ( r x \2 

X^-^w^-u,) = - — ( Po (a) -p (b))- ^ ^ 



3£ ^ 3 " ' , y K " (Vba - W - l/T 2 ) [1/T 2 + (Qba - aOTO/Ii) 



= 8Ne6 (n (a) o (h\) & 

5e ft 5 [P0[ ' m >! (Vba -u- i/T 2 ) [1/T 2 + (Qba - wJTO/T^ ' 

as being the first contributions to the two-level polarization density, including up to fifth order 
interactions. For a summary of the non-resonant forms of the susceptibilities of two-level systems, 
se Butcher and Cotter, Eqs. (6.71)-(6.73). 



4 In the explicit expressions for the linear susceptibility, for example Butcher and Cotter's 
Eqs. (4.58) and (4.111) for the non-resonant and resonant cases, respectively, there are two terms, 
one with Q ba — w in the denominator and the other one with Qb a + uj. The reason why the second 
form does not appear in the expression for the field corrected susceptibility, as derived from the 
Bloch equations, is that we have used the rotating wave approximation in the derivation of the 
final expression. (Recapitulate that in the rotating wave approximation, terms with oscillatory de- 
pendence of cxp[i(Qb a +Lo)t] were neglected.) As a result, all temporally phase-mismatched terms 
are neglected, and in particular only terms with Q ba - w in the denominator will remain. This, 
however, is a most acceptable approximation, especially when it comes to resonant interactions, 
where terms with Q ba -w in the denominator by far will dominate over non-resonant terms. 
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Appendix: Notes on the numerical solution to the Bloch equations 

In their original form, the Bloch equations for a two-level system are given by Eqs. (1) as 

-Av-u/T 2 , 
Au + f3(t)w-v/T 2 , 
-p(t)v-{w-wo)/T x . 
By taking the time in units of the transverse relaxation time T 2 , as 

t = t/T 2 , 

the Bloch equations in this normalized time scale become 

-AT 2 v - u, 

AT 2 u + (3{t)T 2 w - v, 

-0(t)T 2 v - (w - woW/Tl 

In this system of equations, all coefficients are now normalized and physically dimensionless, ex- 
pressed as relevant quotes between relaxation times and products of the Rabi frequency or detuning 
frequency with the transverse relaxation time. Hence, by taking the normalized parameters 

S = AT 2 , 
7(t) - 0(t)T 2 , 
V = Ti/T 2 , 

where S can be considered as the normalized detuning from molecular resonance of the medium, 
-f(t) as the normalized Rabi frequency, and r\ as a parameter which describes the relative impact 
of the longitudinal vs transverse relaxation times, the Bloch equations take the normalized final 
form 

— Sv — u, 
Su + -f(t)w — v, 
-■y(t)v -(w- w )/n. 

This normalized form of the Bloch equations has been used throughout the generation of graphs 
in Figs. 2-11 of this lecture, describing the qualitative impact of different regimes of resonance and 
relaxation. The normalized Bloch equations were in the simulations shown in Figs. 2-11 integrated 
by using the standard routine 0DE45O in MATLAB. 



du 
lit 
dv 
~dt 
dw 
~dt 



du 
dl 
dv 
aW 
dw 



du 

ah 
dv 

dr 
dw 

dr 
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